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Kramers problem for nonequilibrium current-induced chemical reactions 

Alan A. Dzhioe\0 and D. S. Kosoxj^ 

Department of Physics, Universite Libre de Bruxelles, Campus Plaine, CP 231, Blvd. du Triomphe, B-1050 Brussels, 
Belgium 

We discuss the use of tunneling electron current to control and catalyze chemical reactions. Assuming the 
separation of time scales for electronic and nuclear dynamics we employ Langevin equation for a reaction 
coordinate. The Langevin equation contains nonconservative current-induced forces and gives nonequilibrium, 
effective potential energy surface for current-carrying molecular systems. The current-induced forces are 
computed via Keldysh nonequilibrium Green's functions. Once a nonequilibrium, current-depended potential 
energy surface is defined, the chemical reaction is modeled as an escape of a Brownian particle from the 
potential well. We demonstrate that the barrier between the reactant and the product states can be controlled 
by the bias voltage. When the molecule is asymmetrically coupled to the electrodes, the reaction can be 
catalyzed or stopped depending on the polarity of the tunneling current. 



When a molecule is attached to two metal electrodes 
with different chemical potentials or when it is placed on 
a surface under scanning tunneling microscope (STM) 
tip, the electron current fiows through it. It brings the 
molecule out of equilibrium and changes its electronic, 
vibrational and mechanical properties. The interaction 
of nonequilibrium current-carrying electrons with nuclear 
degrees of freedom may catalyze certain chemical reac- 
tions which are not possible under equilibrium condi- 
tions. This opens a possibility to use the tunneling molec- 
ular junction as a nanoscale chemical reactor rather than 
an electronic element of a circuit. The recent experimen- 
tal work has demonstrated that the chemical bonds can 
be selectively broken or formed by the tunneling electron 
current ji"— The tunneling current can even initiate chem- 
ical reactions in the reactants which are brought close to 
each other under the STM tip4.iS. To predict the out- 
come of current-induced chemical reactions and to guide 
the experimental work, we need to develop intuitively 
simple reaction rate theory which takes into account the 
tunneling molecular junction conditions. 

Chemical reactions in a complex environment are tra- 
ditionally modeled as an escape of a Brownian particle 
from a potential well. One usually begins with the (gener- 
alized) Langevin equation for the reaction coordinate and 
then computes the rate at which an effective Brownian 
particle escapes from the potential well over a potential 
barrier (so-called Kramers problem) i^ii However, when 
the molecule is driven out of equilibrium by the tunnel- 
ing flow of electrons through it, the energy or free energy 
surface can not be defined. How can one formulate and 
solve Kramers problem in this case? That is one of the 
key questions which we address in this paper. 

Let us consider a molecule attached to two metal elec- 
trodes. One electrode can be, for example, the metal 
surface and the other one is a STM tip. The molecule is 
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modeled by one electronic spin-degenerate molecular or- 
bital with energy e{x), which depends on some reaction 
coordinate x and the gate voltage. The reaction coordi- 
nate is considered to be a classical variable with corre- 
sponding momentum p and reduced mass m. The nuclear 
Coulomb repulsion energy is U{x). Then the molecular 
Hamiltonian (we use atomic units throughout the paper) 
is 



Hm = e{x) 



P 

2m 



U{x). 



(1) 



Here a\{aa) creates (annihilates) an electron with the 
spin cr in the molecule. The total molecular junction 
Hamiltonian consists of the molecular Hamiltonian ([IJ, 
the Hamiltonians for noninteracting left and right elec- 
trodes, and the molecule-electrode interaction: 



H — Hm + 



a,keL,R 



£kO,l^i~aak 



E 

a,keL,R 



{tka^^aa + h.c). 



(2) 

where ^^/^(acrfe) creates (annihilates) an electron in the 
state ak of either the left (L) or the right [R) elec- 
trodes. Electron creation and annihilation operators sat- 
isfy standard fermionic anticommutation relations. Tun- 
neling coupling matrix element is ■ Since the screening 
length in the metallic electrodes is very short, we assume 
that the voltage bias drops on the interface. Therefore 
there is no external electric field in the molecule. 

To derive Langevin equation for the reaction coordi- 
nate we partition the molecular junction into two parts: 
a "system" and a "bath". The system is the reaction 
coordinate x and the bath consists of all electronic de- 
grees of freedom in the molecule and electrodes. The 
bath degrees of freedom can be projected out from the 
equations of motion and affect the reaction coordinate 
only through effective forces.® Furthermore, we assume 
that the time-scales for electronic and nuclear motions 
can be separated: The electronic degrees of freedom are 
much faster than the motion of the molecule along the 
reaction coordinate x. Therefore, we can assume that 
the electronic steady state is instantaneously established 
along the reaction coordinate trajectory x = x{t). So the 
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nonequilibrium electronic density matrix p{x) depends 
on time only through the parametrical dependence on 
X. The resuh is the Langevin equation for the reaction 
coordinate^"— 

mx = -1:t[p{x)—] -Cx + Sf{t). (3) 

Here the conservative part of the force is given by the 
nonequihbrium analogue of the Hellmann-Feynman theo- 
rem Tr[p{x)^]r^ and Sf{t) (fluctuating force), (i (fric- 
tional force) are nonconservative contribution originated 
from the integrated out electronic degrees of freedomi^ 
The nonconservative forces describe Joule heating, i.e. 
the energy loses due to particle-hole excitations in the 
molecule and metal electrodes. For particular systems, 
the conservative current-induced forces can be com- 
puted with the use of nonequilibrium Green's functions 
within tight-binding approximation or density functional 
theoryi^ii^. 

The noise is taken in the Gaussian form and related to 
the viscosity by the fluctuation-dissipation relation with 
some temperature T: 

(Sfit)) = 0, {Sfit)Sf{t')) = 2CTS{t - t'). (4) 

We assume that the electrons and vibrations have the 
same temperature T although we fully appreciate that 
far from equilibrium vibrational temperature can devi- 
ate from electronic temperature^^"— . Since dH/dx — 
e'{x) a^Oo- -I- U'{x), Eq. ^ becomes 

mx = -e'{x)n{x) - U'{x) ~ (x + 5f{t) ~ (x, (5) 

where n{x) = 2TT[p{x)al.aa] is the nonequilibrium pop- 
ulation of the electronic level in the molecule. 

To complete Langevin equation for the reaction coordi- 
nate ^ we need to know the explicit expression for n(x). 
For a given value of x the population n{x) can be com- 
puted by Keldysh nonequilibrium Green's functions. 
The derivations are relatively straightforward, so we just 
outline them here without giving the full detailsJ^ One 
begins with Keldysh contour-ordered Green's function 
and writes the Dyson equation for it. Applying the Lan- 
greth rules for analytical continuations the Dyson equa- 
tion is solved for and the nonequilibrium population and 
electron current are associated with (lesser) Green's 
function on the real time axis. It results in the following 



expression 
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J TT [uj - s[x) - A{uj)y + {T{uj)y 

Here F = F^+Ffl and A — A^+Aj^ determine the broad- 
ening and shift of the molecular level due to coupling to 
the electrodes. They are given by the real and imaginary 
parts of the electrode self-energy 



l^fcl 



The function /L,i?(w) = [1 + e("-''^'«)/^]"i is the Fermi- 
Dirac electron distribution in the left and right elec- 
trodes. 

The corresponding electron current also depends para- 
metrically on the reaction coordinate and can be readily 
calculated by the Landauer formula.^* In the equilibrium, 
i.e. when the chemical potential of the left electrode 
equals to the chemical potential of the right electrode, 
the current is zero and Langevin equation ([5]) remains 
the same but n{x) becomes the equilibrium electronic 
population noq(a;), which can be computed from © by 
setting /l = Jr. 

Let us introduce the time-dependent probability dis- 
tribution for the reaction coordinate F — F{x,p,t). 
The Langevin equation ^ is equivalent to the standard 
Fokker-Planck equation for the distribution function: 



dt m dx dp '^^^ m dp^ 
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(7) 



Here the effective nonequilibrium potential energy sur- 
face is defined via the integration of the nonequilibrium 
force in Langevin equation ([5]) 

Ueff{x) = U{x)+ I dye'{y)n{y). (9) 

The choice of a;o is not relevant, since the Fokker-Planck 
equation does not depend on it. To elucidate effects re- 
lated to the electron current we separate the effective 
potential into two parts: 

Ueff{x)^U,^{x)+ f dye'{y)Aniy). (10) 

Jxo 

Here the first term describes the equilibrium potential 
energy surface and the second term gives the nonequi- 
librium corrections. The equilibrium potential energy 
surface, Ucq{x), is the potential energy of the molecule 
in the absence of electron current and it includes equi- 
librium charge transfer between the molecule and metal 
electrodes. Since the derivative e'{y) does not depend on 
the applied voltage and on the coupling to the electrodes, 
the nonequilibrium correction is only due to the varia- 
tion of the molecule population An(x) = n{x) — neq{x) 
caused by the current flow. As one sees from Eq.® the 
nonequilibrium correction An{x) is most significant when 
the molecular level is in resonance with the Fermi energy 
of the electrodes, i.e., e{x) w p,L,R- In this case the 
current reaches its maximal value. Moreover, at asym- 
metric molecule-electrode coupling, F^ 7^ F^, An(x) can 
be both positive and negative, depending on the polarity 
of the applied voltage bias. Therefore, the current fiows 
through the molecule can locally increase or decrease the 
nonequilibrium potential energy surface. 

Let us understand the behavior of the effective po- 
tential energy surface based on numerical calculations. 
We take £{x) in the form of bonding orbital for H2 
moleculei^. For simplicity, we assume that both elec- 
trodes have the constant density of states, "qn and ?7l. 
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FIG. 1. Effective nonequilibrium potentials Ueff{x) as functions of the reaction coordinate. The voltage bias is V and the 
asymmetry coefficient g = Tl/^r- 



and they are characterized by the same half-bandwidth 
D = 50. We also assume that the tunneling coupling tk is 
real and independent of k. Then the real and imaginary 
parts of the electrodes self-energy are 



AL,i?(w) 
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L,R 



{lo) 



L,R 



Q{D- 



(11) 

where Tl u = irtj^ r11l,b.- The total width of the molec- 
ular level r = Ti -t-Fij = 0.05 is fixed in our calculations 
but we vary relative contributions of left and right elec- 
trodes. The ratio g = Tl/Tji is the asymmetry coeffi- 
cient. The asymmetry coefficient can be experimentally 
controlled by changing the relative strength of the cou- 
pling tk (HI) of the molecule to left/right electrodes. This 
for example can be accomplished by selecting chemically 
different left and right molecular-electrode Iinkers22i2i or 
simply by changing the distance between the molecule 
and one of the electrodes. The latter can be easily ac- 
complished on STM experiments. 

Fermi energies of both electrodes coincide with the 
band centers, Ef — 0, and the chemical potentials are 
functions of the external applied voltage fiL.R^ = ±0.5^. 
The calculations are performed at the room temperature 
T = 300K. The value of the gate voltage is chosen in 
such a way that e{xmin) = e/j where Xmin = 2.493 deter- 
mines the minimum of H2 ground state energy.— So the 
molecular level is in resonance with the Fermi energy of 
the electrodes. This particular choice of the gate voltage 
is not critical, all our results will qualitatively remains 
the same for other values of the gate voltage. When 
X < Xniin the molecular orbital is below the Fermi energy 
of the electrodes, and when x > x^in it lays above it. 

Fig. [1] shows effective nonequilibrium potential Ueff as 
a function of reaction coordinate x for different values of 
the applied voltage V and the asymmetry coefficient g. 
As one can see from Fig. [TJ the height of the potential 
barrier, AUeff — Ueff{b) — Ueff (a), between product 
and reactant states can be made smaller or larger as the 
current flows through the molecule. For the symmet- 
ric coupling to the left and right electrodes, the barrier 
is always decreased by the voltage bias. This effect is 
mainly caused by the reduction of the electronic popu- 



lation. The asymmetric case is much more interesting. 
Here, the effect of the barrier reduction can be ampli- 
fied or reversed depending on the polarity of the applied 
voltage bias. We consider the case when the molecule is 
coupled stronger to the right electrode {g < 1). When 
the applied voltage is positive V > 0, i.e., the left chem- 
ical potential is larger than the right chemical potential, 
molecular electrons are more depleted than in the sym- 
metric case. Moreover, as follows from Fig.[TJ the applied 
voltage shortens the distance Ax = & — a in reaction coor- 
dinate between the minimum and the top of the barrier of 
the potential energy surface. Both effects collectively re- 
duce the potential barrier height AC/g// stronger than in 
a case of the symmetric coupling to the electrodes. For 
large voltage V = 0.2 and notable asymmetry g = 0.1 
the barrier completely vanishes. It is interesting that by 
reversing the voltage bias we can increase the barrier be- 
tween the product and reactant states. For example, as 
one can see from Fig. [l] when the left electrode is neg- 
atively biased {V = —0.2) and weakly coupled to the 
molecule {g = 0.1) the barrier is increased by 70% as 
compared to the equilibrium. We emphasize that there 
is no electric field across the molecule from the voltage 
bias in our model, therefore the observed physical behav- 
ior is solely due to tunneling current and corresponding 
nonequilibrium changes in molecular electronic popula- 
tion. 

Let us now compute the rates for current induced 
chemical reactions. Since we have already computed 
nonequilibrium effective potential energy surface, we 
can use standard reaction rate theory for our calcula- 
tions. We consider separately two cases: overdamped 



( C ^ Co, where Co — 2yTOC/^q(a) ) and underdamped 

(C ^ Co)- In the overdamped limit, Fokker-Plank equa- 
tion (|5]) becomes one-dimensional Smoluchowski equa- 
tion for the probability distribution and the reaction rate 
k to overcome the barrier can be computed by numerical 
integration^. For small T the reaction rate can be evalu- 
ated analytically by performing the quadratic expansion 
of the effective potential Uef / near the minimum and the 
maximum of the barrier?^ It results into the standard 
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FIG. 2. Nonequilibrium reaction rates computed for vari- 
ous values of the asymmetry coefficient g as a function of 
the apphed voltage. Reaction rates are computed in un- 
derdamped (red curves) C, = O.l^o and overdamped (black 
curves) ^ = lO^o regimes. 



In conclusion, we have demonstrated that the electric 
current which flows through the molecule can be used 
to control chemical reactions. We combined Langevin 
equation for a reaction coordinate with Keldysh Green's 
function calculations of the current induced forces and 
demonstrated how the nonequilibrium, current-depended 
potential energy surface can be defined. The chemical 
reaction is modeled as an escape of a Brownian particle 
from the potential well (Kramers problem). The barrier 
between the reactant and the product states can be con- 
trolled by the bias voltage. We demonstrated that when 
the molecule is asymmetrically coupled to the electrodes 
the reaction can be catalyzed or stopped depending upon 
the direction of the electric current which flows through 
the molecule. 
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expression for the reaction rates 



-AU,ff/T 



(12) 



In the underdamped limit, the reaction rate can be also 
computed analytically^ and it is given by the following 
expression 



CAf/< 



mT 



eff ^-AU,ff/T_ 



(13) 



Using Eqs. (|12ll3p we calculate the reaction rate fc as a 
function of the applied voltage bias for different values of 
g < 1. The results are shown in Fig. [2j For symmetric 
case {g — 1) the reaction rate does not depend on the 
polarity of voltage and it increases with applied voltage. 
The situation changes dramatically when the molecule is 
asymmetrically coupled to the electrodes. The reaction 
is "catalyzed" if the voltage is positive and slows down 
if we reverse the direction of the voltage. It is inter- 
esting that the voltage dependence of the reaction rates 
are almost identical for overdamped and underdamped 
cases. This can be easily understood from the following 
geometrical consideration. Depending on the polarity, 
the tunneling current makes the potential barrier AC/g/ / 
larger or smaller. If the potential barrier increases or de- 
creases, the frequencies U'^^{a) and —U'^jj[b) increase 
or decrease, respectively, too. It leaves the ratio be- 
tween overdamped and underdamped reaction rates al- 
most voltage independent. 
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